function logdet_V = wishart_exp_logdet(W, v)

D = size(W, 1);
args = (v+1-(1:D))/2;
logdet_V = sum(digamma(args)) + D*log(2) + logdet_chol(W);
